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. Abstract 
O. 

^ . We present predictions of the total cross sections for pair production of squarks 

and gluinos at the LHC, including the stop-antistop production process. Our 
calculation supplements full fixed-order NLO predictions with resummation of 
^ ■ threshold logarithms and Coulomb singularities at next-to-leading logarithmic 

(NLL) accuracy, including bound-state effects. The numerical effect of higher- 
order Coulomb terms can be as big or larger than that of soft-gluon corrections. 
For a selection of benchmark points accesible with data from the 2010-2012 LHC 
runs, resummation leads to an enhancement of the total inclusive squark and 
gluino production cross section in the 15-30% range. For individual production 
processes of gluinos, the corrections can be much larger. The theoretical uncer- 
tainty in the prediction of the hard-scattering cross sections is typically reduced 
to the ±10% level. 



1 Introduction 



Despite the good agreement of the Standard Model (SM) with a weahh of experimental 
data, both empirical reasons (e.g. the observation of dark matter) and theoretical argu- 
ments (such as the naturalness problem and the desire for gauge coupling unification) point 
to physics beyond the SM. One of the most thoroughly studied extensions of the SM is Su- 
persymmetry (SUSY). In particular the Minimal Supersymmetric Standard Model (MSSM) 
with i?-parity conservation and superpartner masses at the TeV scale could provide a reso- 
lution to the above issues. The search for SUSY at the TeV scale is therefore a central part 
of the physics program of the Large Hadron Collider (LHC) at CERN. In the context of the 
MSSM, and of any other i?-parity conserving model, the supersymmetric partners of the 
SM particles are produced in pairs, and squarks and gluinos, coupling strongly to quarks 
and gluons, have typically the highest production rates. Experimental searches for SUSY 
have been performed at LEP, the Tevatron and the LHC in various final-state signatures, 
see P for a recent review. For squark and gluino-pair production the tightest bounds 
generically arise from jets-|-missing energy signatures, where the two LHC experiments 
have set lower bounds on the mass of squarks and gluinos of about 800 GeV-1 TeV [2l|3], 
depending on the precise underlying theoretical model assumed. Once upgraded to its 
nominal energy of 14 TeV the LHC should be sensitive to squark and gluino masses of up 
to 3 TeV m. 

SUSY searches and, if squarks and gluinos are discovered, the measurement of their 
properties rely on a precise theoretical understanding of the production mechanism and on 
accurate predictions of the observables used in the analysis. From the theoretical point of 
view, the simplest of such observables is the total production cross section, on which we will 
focus in this work. For QCD mediated processes, such as squark and gluino pair production, 
the Born cross section is notoriously affected by large theoretical uncertainties, such that 
the inclusion of at least the next term in the expansion in the strong coupling constant is 
mandatory for a reliable prediction. Next-to-leading order (NLO) SUSY QCD corrections 
for production of squarks and gluinos were computed in |S] and implemented in the program 
PROSPINO PE]. The corrections are large, up to 100% of the tree-level result, and lead to a 
significant reduction of the scale dependence of the cross section. Electroweak contributions 
were also investigated [HHTT], but found to be much smaller than the QCD contributions, 
less than 5 — 10% of the Born result. 

The size of the 0{as) corrections raises the question of the magnitude of unknown 
higher-order QCD corrections, and makes it desirable to include at least the dominant 
contributions beyond NLO. It is known that a non-negligible part of the full NLO correc- 
tions arises from the partonic threshold region, defined by the limit (3 = ^Jl - AM'^/s 0, 
with M the average mass of the particles produced and s the partonic centre-of-mass en- 
ergy. In the threshold region the partonic cross section is dominated by soft-gluon emission 
off the initial- and final-state coloured particles and by Coulomb interactions of the two 
non-relativistic heavy particles, which give rise to singular terms of the form ag In^'^ /3 and 
as/ (3, respectively. These corrections can be resummed to all orders in as, thus leading to 
improved predictions of the cross section and smaller theoretical uncertainties. Note that 
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to obtain the total hadronic cross section, the partonic cross section is convoluted with 
parton luminosity functions. The convolution scans over regions where (3 is not necessarily 
small, unless M is close to the hadronic centre-of-mass energy s. Hence, in these regions 
the threshold-enhanced terms cannot be expected a priori to give the dominant contribu- 
tion to the cross section. However one often finds that they approximate the exact result 
reasonably well even in regions where /3 is not very small, so that threshold resummation 
is also relevant for improving predictions of the hadronic cross section. 

Resummation of soft logarithms for squark and gluino production at the next-to-leading 
(NLL) logarithmic accuracy have been presented in |T2] - [T6] using the Mellin-space resum- 
mation formalism developed by [I7H20]. Recently the same formalism has been extended 
to NNLL order [21j and applied to squark-antisquark production [22j. These works do 
not resum Coulomb corrections to all orders, though the numerically dominant terms 
are accounted for at fixed order. All-order resummation of Coulomb contributions and 
bound-state effects for squark-gluino and gluino-gluino production were on the other hand 
investigated in [23H25], without the inclusion of soft resummation. In partial NNLL 

resummation of soft logarithms has been used to construct approximated NNLO results 
for the squark-antisquark production cross section. Recently a new formalism for the com- 
bined resummation of soft and Coulomb corrections has been developed [28||29] . and applied 
to NLL resummation of squark-antisquark production ^9\, and NNLL resummation of ti 
hadroproduction [30] . Contrary to the traditional Mellin-space formalism, in our approach, 
which is based on soft-coUinear effective theory (SCET) and potential non-relativistic QCD 
(pNRQCD), resummation is performed directly in momentum space via renormalizat ion- 
group evolution equations [3T1433] . The combined soft-Coulomb effects have been found to 
be sizeable for the case of squark-antisquark production [29] and lead to a reduction of the 
scale uncertainty, as has been observed as well in [22]. 

In this work we extend the results given in [22] to the remaining production processes 
for squarks and gluinos at NLL accuracy, i.e. squark-squark, squark-gluino and gluino- 
gluino production. We also consider separately the production of pairs of stops, which 
requires the extension of the formalism presented in [29j to particles pair-produced in a 
P-wave state. The paper is organized as follows: in Section [2] we give an overview of squark 
and gluino production processes, set up the calculation and briefly review the resummation 
formalism we employ, listing the ingredients needed for NLL resummation. The validity 
of the formalism for P-wave induced processes, which is necessary for resummation of 
the stop-antistop cross section, is established in Appendix \M Numerical results for the 
cross sections are presented in Section [31 including predictions for a representative set 
of the benchmark points proposed in [31] and a comparison to results using the Mellin- 
space formalism [Hl|15]. Finally in Section H] we present our conclusions and outlook. 
Explicit expressions for resummation functions appearing in the NLL cross sections are 
provided in Appendix [Bl while Appedix [C] contains some details on the scales used in 
the momentum-space resummation and on our method to estimate ambiguities in the 
resummation procedure. 
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2 NLL resummation for squark and gluino production 



At hadron colliders, the dominant production channels for squarks q and gluinos g are 
pair-production processes of the form 

N1N2 ss'X, (2.1) 

where A^^i 2 denote the incoming hadrons and s, s' the two sparticles. The total hadronic 
cross sections for the processes (12.11) can be obtained by convoluting short- distance pro- 
duction cross sections (jpp/(s, /i/) for the partonic processes 

pp^ss'X, p,p' e {q,q,g}, (2.2) 

with the parton luminosity functions Lppi{T,fi): 

aN^Ni^sS'xis) = J dr ^ Lpp>{T, fXf)d-pp>{TS, fxj) , (2.3) 

"^0 p,p'=q,q,g 

where Tq = AM^/s, with the average sparticle mass 

M = !^^i±^. (2.4) 
The parton luminosity functions are defined from the parton density functions (PDFs) as 

Lpp>{T,IJ,)= dXidx2S{xiX2 - T)fp/Ni{xi,IJ,)fpf/N2{x2,IJ') ■ (2.5) 
^0 

We perform a NLL resummation of threshold logarithms and Coulomb corrections to 
the partonic cross section, counting both a^/ /3 and In /3 as quantities of order one, where 
/3 = (1 — 4M^/s)^/^ is the heavy-particle velocity. Our predictions include all corrections 
to the Born cross section of the schematic form 

^S,^^ oc Yl [j) [ln/3^o(«.ln/3) +^?i(«.ln/3)]. (2.6) 

A resummation at NNLL accuracy in the counting ln/3 ~ 1, a^/ /3 ~ 1 , which is beyond 
the scope of this paper but has recently been performed for top-pair production [30], would 
include in addition terms of the form exp[as5f2(as lu/S)] and corrections of order 0{as) 
relative to the NLL cross section, including NLO corrections to the Coulomb potential and 
other higher-order potentials, as well as the non-logarithmic one-loop hard corrections. The 
recent NNLL calculation of squark-antisquark production [22] included the corrections of 
the (72-type related to soft corrections and the hard 0{as) corrections, but kept only the 
(qj^Z/S) ^-term in the sum over k. In Section 1271] we collect some facts about the production 
processes of gluinos and the superpartners of the light quarks at LO and NLO while the 
formalism employed for the NLL resummation is reviewed in Section 12.21 The production 
of stop pairs is included in 12. 3^ while details about the choice of the soft scale in the 
momentum-space resummation formalism and our procedure to estimate the remaining 
theoretical uncertainty are discussed in 12.41 
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Figure 1: Ratio of the LO production cross sections for the processes (12.71) to the total 
Born production rate of coloured sparticles, asusY, for the LHC with y/s = 7 TeV. Left: 
Mass dependence for a fixed mass ratio rriq = rrig = M. Right: Dependence on the ratio 
mg/rriq for a fixed average mass (m^ + nig)/ 2 = 1.2 TeV. 



2.1 LO and NLO results 

At leading order [35H37], the following partonic channels contribute to the production of 
light-flavour squarks and gluinos: 

9Qi M, 9Qi -> M , 
99, mi 99 , (2.7) 

where i, j = u, d, s, c, b. At NLO further partonic processes contribute to the cross 
section. To keep the notation as simple as possible, in (12. 7p we have suppressed the 
helicity and flavour indices of the squarks. It is understood that in the predictions 
for the cross sections presented below the contributions of the ten light-flavour squarks 
{ul/r, di/R, cl/r, sl/r, bL/R) are always summed over. Furthermore, the ten scalars are 
assumed to be degenerate in mass, with the common light-flavour squark mass given by m,g. 
In the following, the charge-conjugate subprocesses for squark-squark and gluino-squark 
productions will be included in our results. As input for the convolution (12. 3p we will 
use the MSTW08 set of PDFs |3BJ at the appropriate perturbative order (the LO PDFs 
for Born-level predictions and the NLO PDFs for the NLO and NLL results) and set the 
factorization scale to the average mass of the produced sparticles, fif = M. We use a set 
of PDFs with an improved accuracy at large x provided to us by the MSTW collaboration 
that has also been employed for the NLL results in [T6] . 

To illustrate the relative magnitude of the various processes depending on the squark 
and gluino masses, the ratio of the total hadronic cross section for the processes (12. 7p to 
the total inclusive cross section for squark and gluino production ctsusy = (Tpp^qq+qq+gq+gg 
is shown in Figure [T] for the LHC with ^/i = 7 TeV centre-of-mass energy. From the left- 
hand side plot, showing the relative contributions of the various processes as a function of 
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Figure 2: NLO i^-factor for the processes ( 12. 7p at the LHC with ^/s = 7 TeV. Left: Mass 
dependence for a fixed mass ratio rriq = rrig = M. Right: Dependence on the ratio mg/irtq 
for a fixed average mass {niq + rrig) /2 = 1.2 TeV. 



a common squark and gluino mass, it can be seen that squark-squark and squark-gluino 
production are by far the dominant channels over the full mass range considered. In the 
right-hand side plot, the relative contributions are shown as a function of the squark- 
gluino mass ratio for average mass = 1.2 TeV and it is seen that only for gluinos 
that are significantly lighter than squarks, gluino-pair production becomes the dominant 
production channel. In Figure |2] we show the i^-factor i^NLO = o"nlo / <^lo for the SUSY- 
QCD corrections for the various production processes as obtained from PROSPINO [6] 13 
The corrections are positive and enhance the cross section from 40% for squark-gluino 
production with light sparticle masses up to 100% or larger for squark-antisquark and 
gluino-pair production at large sparticle masses. 

Since the focus of this work is on higher-order corrections that are enhanced in the 
threshold limit /3 — )■ 0, we consider here the corresponding terms appearing at NLO. To 
this end, we decompose the partonic cross section appi in a complete colour basis, and 
parametrize the higher-order corrections as 



CTpp' 



Rat 



The sum is over the irreducible colour representations appearing in the decomposition 
R<^R' = Ra, where R, R' are the SU{?>) representations of the two final-state sparticles. 
The relevant decompositions for squark and gluino production are given by 



(2.9) 

To ©27. 



qq : 


3$ 


53 


= 1 6 


D8, 


qq : 


3$ 


53 


= 36 


D6, 


■ 


3$ 


58 


= 36 


D 6 © 15 , 


99 ■ 


8$ 


58 


= 1 6 


D8, ©8a©10 



^To see the genuine size of the NLO corrections, the i^-factors have been computed using the NLO 
PDFs for the Born cross sections. 
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The explicit basis tensors for the various representations have been constructed in [28] (see 
also [T3l[Tl]), where it has been shown that an s-channel colour basis based on the decom- 
positions (I2.9P is advantageous for the all-order summation of soft-gluon corrections. In 
Eq. ( 12. Sp cTpp/'^" represents the tree-level cross section for a given process in colour channel 

Ra, while the f^J' " are colour-specific NLO scaling functions. The colour-separated Born 
cross sections for squark and gluino production are available in [TBHT^ . The NLO scaling 
functions, on the contrary, are only known numerically in their colour-averaged form [S]. 
However, a simple formula is available for the threshold limit of the NLO scaling functions, 
containing all the threshold-enhanced contributions, for arbitrary colour representation 

Ra in]: 



In' I 1 + 



H J 24 

-A{Cr^ + A{Cr + aO) In (^^) + 12Cij. + + (2-10) 

with M the average mass of the two particles produced (12.41) . while denotes the reduced 
mass, rrir = msTngi / {nis + mg')- C*r, Cr' and Cr^ are the Casimir invariants for the colour 
representations of the initial-state particles, p and p', and for the irreducible representation 
Ra of the SUSY pair. The coefficients Dr^ of the Coulomb potential for the production of 
heavy particles in S't/(3) representations R and R' in the colour channel Ra are given in 
terms of the quadratic Casimir operators for the various representations: 

Dr^ = \{Cr^-Cr-Cr.), (2.11) 

where negative values correspond to an attractive Coulomb potential, positive values to a 
repulsive one. The numerical values for the representations relevant for squark and gluino 
production can be found in [29lll0] and are collected in Table [H The coefficient ^p^/ " is 
the one-loop contribution to the hard matching coefficient appearing in Eq. (I2.13P below, 
and represents the only process-specific quantity in Eq. (12.101) . It has been obtained 
recently for squark-antisquark production and gluino-gluino production |22l|2l] but is not 

fl) R 

known yet for the remaining production processes. The knowledge of h^^i" " is required 

for NNLL resummation [22], but not at NLL accuracy as considered here, so the /ipp/'^" 
will be always set to zero in the following. Using the Born cross sections for the different 
colour channels [ISJIH] and (12.101) one can reproduce the threshold expansions of the NLO 
corrections in [5]o 



^In [5] there is a typo in the sign of the Coulomb correction for like-flavour qq production, i.e. the 
function fqq^^ in eq. (54). Also note that [5] choose to expand the cross section of qg production in the 
variable /3 — -^/l — Amqirig /{s — {rrig — nrigY) « (3^J M /2mr which leads to the appearance of additional 
\n{mr / M) terms. For this process they also observe an apparent non-factorization of the colour-averaged 
NLO threshold corrections from the Born cross section. This is nevertheless consistent with (j2.10|) since the 
Born qg cross section in the colour-triplet channel is not proportional to that for the other channels [14) . 
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Table 1: Numerical values of the coefficients of the Coulomb potential (12. lip for squark 
and gluino production processes. Negative values correspond to an attractive potential. 

In Figure [3] we study to which extent the full NLO corrections as obtained from 
PROSPINO are approximated by the singular NLO corrections, obtained by dropping all 
constant terms from fl2.10p . including In 2 terms, and convoluting the resulting partonic 
cross section (12. 8p with the parton luminosities. For the Born cross sections a^J,'^" in (12. 8p 
the exact expressions, without use of the threshold approximation, have been kept, but 
colour channels with a vanishing threshold limit of the Born cross section at leading or- 
der in /3 have been dropped. For the case of degenerate squark and gluino masses it is 
seen that the difference of the threshold-enhanced contributions to the full NLO correc- 
tions is at the 10 — 30% level over the whole mass range considered, with the exception of 
the squark-squark production channel where the threshold contributions account for only 
40 — 50% of the full NLO corrections. For rrig > rriq the singular terms overestimate the 
corrections for the processes involving gluinos, while the agreement for squark- ant isquark 
production improves. For rrig < rriq the singular terms approximate the full corrections 
very well for all processes apart from squark-squark production. Comparing to Figured], it 
is seen that the singular terms capture the NLO corrections to the dominant processes for 
larger mass ratios (i.e. squark-squark production for rrig ~ 2mq and gluino-pair production 
for rrig ~ O.brriq) rather well. For degenerate squark and gluino masses the quality of the 
threshold approximation for the dominant squark-squark and squark-gluino processes is 
somewhat worse. In all cases, the inclusion of the threshold enhanced NLO corrections in 
addition to the Born terms improves the agreement with the full NLO results. This moti- 
vates the computation of the higher-order threshold-enhanced terms through resummation, 
as performed in the remainder of this work. 

2.2 Soft-gluon and Coulomb resummation 

Next, we briefly review the formalism for the combined resummation of soft- and Coulomb- 
gluon corrections [SHUSH] and provide the relevant ingredients for squark and gluino pro- 
duction at NLL accuracy. We also discuss some features of our implementation that differ 
from that used previously for squark- antisquark production in [29] . 

The combined soft-Coulomb resummation for the production of squarks and gluinos is 
based on a factorization of the hard-scattering total cross section for partonic subprocesses 
of the type (12. 2p . It can be shown that near the partonic threshold, 

s {rris + rris^y , (2.12) 

the partonic cross-section factorizes into three contributions [29], a hard function H, a 



7 



2.0 - 
1.5 - 

1.0 ; 

0.5 = 



m - =m - 

S 9 



11 

99 

qg 



Ao" NLOsing/Ao" NLO 

2.0 r 



600 800 1000 1200 1400 1600 1800 2000 



M[GeV] 



1.5 



1.0 



0.5 



q'q 

■ gg 



= 1.2 TeV 



1.0 



1.5 



2.0 



Figure 3: Ratio of the singular NLO contributions obtained from (12.101) to the exact NLO 
corrections for the LHC with ^/s = 7 TeV. Left: Mass dependence for a fixed mass- 
ratio rriq = rrig = M. Right: Dependence on the ratio rrig/mq for a fixed average mass 
{mq + mg)/2 = 1.2 TeV. 



soft function W containing soft gluons to all orders, and a potential function J summing 
Coulomb gluon exchange: 

V(5,/^) = $^i^p^?K-,^5,/i) J duJR^{E-'^)W''-{u,fi). (2.13) 

Ra 

Here E = \/I—2M is the energy relative to the production threshold and the sum is over the 
colour representations (12. 9p . In (I2.13P the s-channel colour basis mentioned above, that can 
be shown to diagonalize the soft function to all orders |28], is chosen for the hard-scattering 
amplitudes. Independent of the sparticle type, the soft function then depends only on the 
colour representations of the initial-state partons and the irreducible representation of 
the sparticle pair appearing in the decompositions (12. 9p . in agreement with the picture 
that soft-gluon radiation is only sensitive to the total colour charge of the slowly moving 
sparticle pair [20]. The formula (I2.13P has been derived in [2^ for S-wave dominated 
production processes up to NNLL accuracy. This covers all production processes of squarks 
and gluinos, apart from quark-antiquark initiated stop-antistop production, that proceeds 
through a P-wave. The applicability of the formalism to stop production is discussed in 
Section 12.31 and Appendix [XI 

It can be argued that the natural scale for the evaluation of the hard function in (I2.13p . 
leading to well-behaved higher-order corrections, is of the order of fih ~ 2M, while the 
natural scale for soft-gluon radiation is of the order of fig ~ We use the momentum- 

space resummation formalism of [311433] to evolve the soft and hard functions from their 
natural scales to the factorization scale fif used for the evaluation of the parton distribution 
functions, commonly taken to be of the order /i/ ~ M. In this way, logarithms of fis/fJ'f ~ 

are summed to all orders. The precise prescription for the choice of the soft scale 
adopted in our calculation is discussed in 12.41 The exchange of multiple Coulomb gluons 
can be summed up using the method of Coulomb Green's functions in non-relativistic 
QCD [1M3]. 
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For resummation at NLL accuracy, the leading-order hard and soft functions are re- 
quired as fixed-order input to the evolution equations. The leading-order soft function 
is trivial, W^^^^°'{u) = 6{u). The leading-order hard functions are obtained from the 
threshold limit of the Born cross section for the colour channel Ra [29] , 



Although the Born-cross sections in the threshold limit appear on the left-hand side 
in fl2.14p . we keep the exact expressions in our numerical implementation, so the hard func- 
tions in practice depend on s. This incorporates some higher-order terms in P, albeit not 
systematically. |fl Here our current treatment differs from that used for squark-antisquark 
production in [29] where only the threshold limit of the Born cross section was used to 
compute the hard function. 

For the resummation of Coulomb corrections, we use results for the non-relativistic 
Coulomb Green's function obtained for top-quark production at electron-positron colliders 
and stop production [41j,|44j. For positive values of E and vanishing decay widths of the 
sparticles, the S-wave potential function is given by the Sommerfeld factor 

JrAE) = ^ \^ [e^^n^-s^— _ ij ^ ^>o, (2.15) 

with the coefficients of the Coulomb potential given in 02.111) . For an attractive 
potential, a series of bound states develops below threshold with energies 

En = (2.16) 

Their contribution to the S-wave potential function is given by 

Jl°^-\E) = 2f2KE-Er,) ^ 2m,(-D^Ja, y ^ 17) 

n=l ^ ' 

For sufficiently broad squarks and gluinos with decay widths exceeding the binding energy 
of the would-be bound states, Vg + Vgi > \Ei — -E2I, the bound state poles are smeared out 
by the finite lifetimes0 We consider here a situation where the widths of the squarks and 
gluinos are large enough to prevent the formation of bound states, but small enough that 



■^This corresponds to the treatment of [TMIS] . up to the fact that we set the hard function for a given 
production and colour channel to zero if the Born cross section vanishes at threshold at leading order in ^, 
even if the full Born cross section for this channel is non- vanishing. This affects the subprocesses qq — gg 
in the singlet channel and qiqi — > qiqi in the triplet channel. The numerical effect is however negligible. 

''The actual formation of bound states is possible if the decay widths of the sparticles are smaller than 
the decay-rate of the bound state. For long-lived gluinos, or light stops without allowed two-body decays, 
higher-order corrections to bound-state production have been recently obtained in |45p46j . Since the bound 
states decay into gluon or photon pairs, this scenario leads to very different collider signatures compared 
to the missing-energy signatures of continuum production and is not considered further here. 
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the use of a narrow- width approximation is justified, which is the case for SUSY scenarios 
with moderate mass ratios of squarks and gluinos where Ts/Mg ~ 1%. The contributions 
to the total cross section below the nominal production threshold, s = 2M, can then be 
included by setting the sparticle widths to zero and including the bound-state poles (12 .17^ . 
For other observables, the finite width can be taken into account to a first approximation by 
the replacement E — )• E + i(Ts + Tsi)/2 in the potential function, see e.g. j23]425] for recent 
studies of the invariant-mass spectrum of gluino-pair and squark-gluino production. The 
study of finite-width corrections for larger decay widths (e.g. for gluino masses rrig > 2mq) 
is left for future work. In the numerical results presented in this work, the contributions 
of the bound state poles for E < will always be included in our default implementation, 
and are convoluted with the resummed soft function as described in [30]. Note that in the 
previous results for squark-antisquark production [29] the bound-state corrections have 
been added without soft-gluon resummation. 

The resummed cross section at NLL accuracy is obtained by inserting the potential 
function (I2.15P and the solutions to the evolution equations of the hard and soft func- 
tions [28|H7] into the factorization formula (I2.13p . Using the solutions in momentum-space 
obtained in |33], the NLL cross section is written as 



(2.18) 

Here the label i jointly refers to the colour of the initial-state partons and the represen- 
tation Ra of the sparticle pair. The function rj = '^""C^^+'^r') + ... contains 
single logarithms, while the resummation function Ui sums the Sudakov double logarithms 
a^log^ and a^log^ Explicit expressions up to NLL accuracy are given in Appendix 

IBI For fig < fif the function rj is negative and the factor u'^'^~^ in the resummed cross sec- 
tion (I2.18P has to be understood in the distributional sense, as discussed in detail in [30] . 
We have used the non-relativistic expression E = in the argument of the potential 
function, that is valid near the partonic threshold ( 12.12^ . This follows the default treat- 
ment of top-pair production in [30] and leads to the customary expansion of the cross 
section in /3 (see Eq. I2.10p . We also perform this replacement in the definition of the 
hard functions (12.140 . The difference between this default implementation and the results 
obtained by consistently keeping the expression E = \f§ — 2M (as in the previous results 
for squark-antisquark production [29]) will be used to estimate the effect of subleading 
terms in the cross section, as discussed in Section 12. 4[ 

In order to assess the importance of the Coulomb corrections and to compare to the re- 
sults of the Mellin approach [13] we will also present results without Coulomb resummation, 
obtained by inserting the trivial potential function J^'^^E) = ^"^"^^ \/E/2mr into (12.181) . 
In this approximation, that we will denote by NLL^.!./!, a fully analytical expression for the 
resummed cross section can be obtained: 



NLL.,, ^ ^ aS!(5,/x.)t/.(M,/.„/.„/.,)^^^-^/3^'^. (2.19) 
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Figure 4: Tree-level diagram topologies contributing to qkQi — > qiqj. 

Since contributions to the cross section from outside the threshold region can be numer- 
ically non-negligible, we match the NLL resummed cross section to the fixed-order NLO 
calculation by subtracting the NLO expansion of the NLL expression and adding back the 
full NLO corrections: 

+ (2-20) 

where d'^^^is) is the fixed-order NLO cross section obtained in standard perturbation 

theory, as implemented in PROSPINO [6], and o^^^^^'' is the resummed cross section expanded 
to NLO, as given in [29]. The total hadronic cross section at NLL is then obtained by 
convoluting (12.201) with the parton luminosity, as in (12. 3p . 

2.3 Stop-antistop production 

Beside the channels listed in (12. 7p . in Section [3] we will also present predictions for stop-pair 
production: 

gg^QiQi ijij , (2.21) 

where only the initial states appearing at leading order have been shown. The NLO SUSY- 
QCD corrections have been computed in [7] and are implemented in PROSPINO [6j. Contrary 
to the light-fiavour squark case, in most scenarios the mixing of the two weak eigenstates ti, 
tfi, and the mass difference of the resulting mass eigenstates, h = cos 6'j + tn sin 9^, t2 = 
—ti smOi + tRCOsOi is non-negligible. Off-diagonal production of the mass eigenstates, e.g. 
tit2, appears at NLO in SUSY-QCD, and through electroweak contributions. It is therefore 
suppressed compared to diagonal production |7,48| and will not be considered here. It must 
also be mentioned that because of the absence of a significant top-quark component inside 
the nucleon the processes tt and gt first contribute to the cross section at NLO, and are thus 
numerically suppressed. The NLO i^-factor for the process PP — )■ titi is shown in Figure 
[5] for a centre-of-mass energy of 7 TeV, and the mass range M = 100 — 1000 GeV. As can 
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Figure 5: Left: NLO fT-factor for stop-pair production at ^/i = 7TeV as a function of the 
stop mass. Right: ratio of the singular NLO contributions obtained from Eqs. fl2.10p and 
fl2:25|l to the full NLO cross section for PP iiii. 



be seen, NLO corrections are in the 50 — 60% range. The predictions for the second mass 
eigenstate differ only in the fixed-order NLO results, and for a given mass the numerical 
difference between the cross sections for iiii and ^2^2 production is below 2% for the mass 
range considered in this work. We therefore omit results for the process PP — )■ ^2^2- 

Contrary to the production of a light squark-antisquark pair, stop-pair production in 
the qq channel cannot be mediated by a t-channel diagram like in Figure Hl^a), again due 
to the extreme suppression of top-quark PDFs inside the proton. As a result, at LO in 
QCD a stop-antistop pair is produced in a P-wave state in quark- ant iquark collisions. As 
shown in Appendix \^ the resummation formalism can be extended in a straightforward 
way to gg-initiated stop-antistop production at NLL, and the only modification of f l2.13p 
is the replacement of the potential function by that appropriate for P-wave processes. For 
stable particles, above threshold the result is given by [44J (see also [49j) 

Jl {E) = 2mrE i^l + ^^^1 ^ j Jn^ {E) , E> 0. (2.22) 

The bound-state contributions of the P-wave Green's function can be found in |33] but are 
not needed here, since only the repulsive colour-octet channel appears in our application to 
stop-pair production. By expanding (12. 22^ in the strong coupling constant one obtains the 
coefficients of the fixed-order Coulomb corrections. While the one-loop Coulomb correction 
agrees with that for S-wave production (12.10p . the second Coulomb correction ^^^(vr^-l- 
3) / differs from the S-wave case. Due to the different normalization of the P-wave Green's 
function, the definition of the leading-order hard functions for P-wave production reads 

^(o),i^..^^ _ (2m,)^ // E V (o),r„ , . 

In addition to the the combined soft / Coulomb resummation, we will again consider the 
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NLLs+ft approximation where the trivial P-wave potential function J^°\E) = ^^^(E/2m,)3/2 
is used in the resummation formula, leading to the analytical result 

Cp" = E f^MiM, /i., /i., /^/) 4^^'|T) l^''- (2.24) 

Flty. ^ 

In analogy to the S-wave result fl2.10p . one can use the resummation formalism to 
obtain the threshold-enhanced one-loop scaling functions for P-wave production in the 
colour channel Ra from initial-state partons in the representations r and r': 



fpp',p = — ^vir+^^^^+^'^ 

Rc, + "3" (^ 
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In agreement with [15] one finds that the coefficient of the single logarithms related to 
initial-state radiation is multiplied by a factor of | compared to the S'-wave case while 
the double logarithm and the logarithms related to final state radiation proportional to 
Cr^ are unchanged. In addition, the constant terms are different which is irrelevant at 
NLL accuracy but has to be taken into account if one aims to extract the one-loop hard 
function from a computation of the NLO cross section. In Figure [5] we study the accuracy 
of the threshold approximation defined by inserting the NLO singular terms, obtained from 
Eqs. fl2.10p and f l2.25p by dropping constant terms, into (12. 8p . The ratio of the singular 
NLO corrections to the full NLO corrections to the hadronic cross section obtained using 
PROSPINO is shown in Figure O (right plot). Analogously to squark-antisquark production, 
the threshold terms provide an excellent approximation of the full NLO result. 



2.4 Scale choices 

As explained in Section 12. 2[ the resummed partonic cross section (I2.13P depends on a 
number of scales related to the factorization of hard, soft and Coulomb effects. The 
dependence on these scales would cancel in the exact result, but a residual dependence 
remains at a given logarithmic order. As already pointed out, our default choice for the 
factorization and hard scales are fif = M and fih = 2M, respectively. On the other hand, 
the resummation of all NLL effects related to Coulomb exchange requires that the scale in 
the potential function Jr^ is chosen of the order of \/2m^M [i, which is the typical virtuality 
of Coulomb gluons. A more detailed analysis shows in fact that for an attractive Coulomb 
potential the Coulomb scale freezes when (3 ~ \DR^\as, due to bound-state formation. We 
thus choose the scale in Jr^ to be 

/ic = Max|2a,(/ic)m,|D«J,2v/2m,M/3} . (2.26) 

Note that, for a repulsive potential, D^i^ > 0, no bound states arise, so that fl2.26p is not 
completely justified. However in this case resummation of Coulomb corrections leads to 
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small effects, and Jr^ vanishes for small /3, so that the precise choice of nc in this limit 
has a negligible numerical impact on predictions of the cross section. 

The choice of the soft scale fig presents some subtleties. The exponentiation of all 
NLL ln/3-terms in the partonic cross section would require a choice Hs ~ However 
a running scale leads to strong oscillations of the cross section for small /3, due to the 
prefactor e~'^'^^^ /r{2ri) in fl2.18p . amplified by the factor u"^^ and terms in the function Ui, 
and eventually hits the Landau pole of the strong coupling constant when /3 — )■ 0. To 
overcome these problems two different approaches have been used in the literature: 



Fixed fig'. In pTH55] the choice of a fixed soft scale was advocated. Such a scale is 
determined from the minimization of the one-loop soft corrections to the hadronic cross 
section, 

= /i,— > / dr Lpp^T, fig) — - — — . (2.27) 



In this approach threshold logarithms are resummed in an average sense and not locally 
at the level of the partonic cross section. However one can argue that, for threshold 
dominated processes, the choice fl2.27p preserves the hierarchy between the soft and short- 
distance scales and that logarithmic corrections log(l — tq) to the hadronic cross section 
are correctly resummed. This was the method adopted in [29] for resummation of the 
squark-antisquark production cross section. The explicit value of the scales determined 
with the minimization procedure fl2.27p are given in Eq. (IC.ip in Appendix O 



Running fig'. For the NNLL resummation of soft effects in tt production presented in [30] 
a different approach was adopted. There a running soft scale, 

/i> = , (2.28) 

was used in the interval /3 > (3cnt, and replaced by a fixed soft scale 

/i< = kgM^l, (2.29) 

below the cutoff. With this scale choice, logarithms of /5 are exponentiated locally in the 
partonic cross section in the large-/9 region, where the use of a fixed soft scale cannot be a 
priori justified. On the other hand if Pcut is not too big, in the lower interval the hadronic 
cross section is in fact dominated by logarithms of (3cut, as can be explicitly checked by 
convoluting the partonic cross section with toy parton luminosities [33], so that the use of 
a fixed scale once again correctly resums the dominant logarithms. The precise value of 
Pent is chosen through the prescription described in |30|, which is reviewed in Appendix 
O The default choice for the prefactor kg adopted here is fc^ = 1. We have observed that, 
for the SUSY processes considered here, the NLL expression and its NNLO expansion are 
generally stable against variations of kg for this choicejfl 



^Our choice of ks deviates from the one adopted in [30 , where kg = 2. This corresponds to resumniing 
some of the In 2 terms in the fixed-order cross section alongside the threshold-enhanced In /? contributions. 
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The two possible choices of the soft scale fig just discussed are one of the ambiguities 
associated with threshold resummation. Others are related to the choice of the hard and 
Coulomb scales and to power-suppressed terms which are not controlled by resummation. 
Additionally, one has to consider the ambiguity arising from the choice of the factorization 
scale /i/. The latter clearly also applies to the fixed-order NLO result. Thus, to reliably 
ascertain the residual uncertainty of the fixed-order and resummed results we present in 
Section [3|, we adopt the following procedure: 

• Scale uncertainty: for both the NLO and NLL result the factorization scale is 
varied between half and twice the default value, i.e. M/2 < fif < 2M. For the NLL 
result, this is done keeping the other scales fih, fic A^s and the parameters /3cut and 
kg fixed. 

• Resummation uncertainty: both hard and Coulomb scales are varied between 
half and twice the default values, i.e. M < ji^ < AM and jji'c 1'^ < l^c < '^f^c^ ^ 
where /i^^ is the solution of the implicit equation (12.261) . In addition, for the NLL 
implementation with a fixed soft scale, /i^ is varied between half and twice its de- 
fault value, while for the running-scale implementation uncertainties related to the 
choice of /3cut and kg are estimated according to the procedure given in [30j (and 
reviewed in Appendix [C]). Finally, as anticipated below (I2.18p . we take the difference 
in parametrizing the resummed cross section in terms of /3 or ii^ as a measure of 
the effect of power-suppressed terms. All the scales and the parameters /3cut and kg 
are varied one at the time keeping the other fixed to their central values, and the 
resulting errors are summed in quadrature. 

• PDF uncertainty: we estimate the error due to uncertainties in the PDFs using 
the 68% confidence level eigenvector set of the MSTW08NLO PDFs [38]. 

An additional source of error arises from the uncertainty on the a^-determination. This 
effect has been found to be of the order of 3% for the NLO cross sections of squark- 
squark, squark-antisquark and squark gluino production and up to 8% for gluino-pair 
production [16]. We expect a similar uncertainty of the NLL results. 

In the following we will often refer to the sum in quadrature of scale and resumma- 
tion uncertainty as "total theoretical uncertainty". Note that the terminology adopted 
here differs slightly from the one used for ti production in [30] where the errors from 
variation of the hard and Coulomb scales, and of the soft scale for the fixed-scale imple- 
mentation, had been incorporated into the scale uncertainty, while we consider them as 
resummation ambiguities. Additionally, in [30] independent and simultaneous variations 
of the factorization and renormalization scale have been considered, whereas in this work 
we identify the factorization and renormalization scales and vary them as one scale, i.e. 
M/2 < fif = fir < 2M. This is the default procedure implemented in the numerical code 
PROSPINO used for the computation of the fixed-order NLO result [6]. 

It is interesting to study how the choice of a fixed or running soft scale affects the 
NLL resummed cross section, especially in view of the uncertainties just discussed. In 
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Figure 6: Resummation uncertainty for the NLL resummed result with a running soft scale 
(NLL, solid blue) and a fixed soft scale (NLLfixed, dashed red) for squark-antisquark (top- 
left), squark-squark (top-right), squark-gluino (centre- left), gluino-gluino (centre-right), 
stop-antistop (bottom-left) production and the inclusive gluino and light-flavour squark 
cross section (bottom-right) at LHC with ^/s = 7 TeV. The central line represents the 
i^'-factor for the default scale choice, while the band gives the resummation uncertainty 
associated with the result. See text for explanation. 
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Figure [6] we plot the NLL i^-factor, defined in Eq. (13. 2p . as a function of a common SUSY 
mass M = irig = rrig for the four processes hsted in (12.71) . for a centre-of-mass energy of 
7TeV (the situation for y/s = 14 TeV is quahtatively similar). Results for the stop-pair 
production process and the total SUSY production rate are also shown. The thick lines 
represent the central values for the two implementations, whereas the bands (delimited by 
thinner lines) correspond to the resummation uncertainty as defined above. The central 
values are in good agreement for squark-antisquark and gluino-pair production, and for 
squark-gluino production at larger masses. For squark-squark production the agreement is 
less satisfactory, especially for smaller masses. This is consistent with the observation from 
Figure [3] that the NLO corrections for squark-squark production are not as dominated 
by the threshold contributions as those for the other processes. In all cases, however, 
the two different NLL predictions are consistent with each other once the uncertainty 
associated with the resummation procedure is taken into account. It can also be seen that 
the uncertainty band for the fixed-scale implementation NLLgxed is mostly contained inside 
the uncertainty band of the running-scale result, with the possible exception of the small- 
mass region. In light of this, in Section |3] we will take the matched NLO/NLL result, Eq. 
(I2.20p . with a running soft scale, Eqs. (I2.28P and (12. 29 p . as our default and best prediction. 

3 Numerical results 

In this section we present numerical results for the cross sections of the five SUSY processes 
introduced in Section [2?T] and [2731 In Section [3?T] we discuss the impact of the NLL soft and 
Coulomb corrections on the central value of the total cross sections and the uncertainties for 
the production of light-flavour squarks and gluinos. In Section 13.21 we provide predictions 
for a selection of the benchmark points defined in [M]- The results for stop-antistop 
production are presented in 13.31 In order to facilitate the use of our results, the arXiv 
submission of this paper includes grids with predictions for the LHC with i/s = 7 TeV, 
for light-fiavour squark and gluino masses from 200 — 2000 GeV and stop-masses from 
100 — 1000 GeV. We also provide a Mathematica file containing interpolations of the cross 
sections with ~ 1% accuracy. An update of the grids will be provided in case a higher 
centre-of-mass energy will be available in the 2012 run. 

3.1 Squark and gluino production at NLL 

To illustrate how different classes of corrections contribute to the total cross section, we 
introduce three different NLL implementations: 

• NLL: our default implementation. Contains the full combined soft and Coulomb 
resummation, Eq. (I2.18p . including bound-state contributions below threshold, Eq. 
(I2.17p . For the soft scale we adopt the running scale given in Eqs. (12. 28 p . (I2.29p . 

• NLLno Bs: as above, but without the inclusion of bound-state effects. 
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• NLLs+/i: this implementation includes resummation of soft and hard logarithms only, 
without Coulomb resummation. This is obtained using Eqs. (I2.19P and fl2.24p . 

The three NLL approximations defined above are always matched to the exact NLO results 
computed with PROSPINO, according to (12.201) . As input for the convolution with the 
parton luminosity functions, Eq. O, we adopt the MSTW08NLO PDF set [38] and 
the associated strong coupling constant as{Mz) = 0.1202. Unless otherwise specified, the 
parameter r, defined as 

r = — , (3.1) 

rrig 

is set to one. 

We start presenting results for the NLL i^-factor, defined as 

^matched 

-^NLL = MT,n ' (3-2) 



a 



where a'^^^'^^'^'^ is our matched result for one of the NLL implementations defined in the 
beginning of this Section and cr^^° the fixed-order NLO result obtained using PROSPINO. 
The NLL-i^'-factor for LHC with 7 TeV centre-of-mass energy is plotted in Figure U\ for 
the four light-squark/gluino production processes and the mass range rriq = nig = 500- 
2000 GeV. The results for a/s = 14 TeV and the mass range ruq = rrig = 500-3000 GeV are 
given in Figure E The NLL corrections for our default implementation (solid blue lines) 
can be large, with corrections to the fixed-order NLO results of up to 120% in the upper 
mass range for gluino-gluino production at 7 TeV. The higher-order effects are smaller, but 
still sizeable, for the other three processes, due to the smaller colour charges involved in 
squark-antisquark, squark-squark and squark-gluino production. Furthermore, for a fixed 
SUSY mass the i^'NLL-factor decreases from 7 to 14 TeV, consistently with the expectation 
that at lower centre-of-mass energies the threshold region plays a more prominent role. 

The effect of including Coulomb resummation and its interference with soft resumma- 
tion is on average as large as (or even larger than) the effect of pure soft and hard cor- 
rections, as can be seen comparing our default implementation NLL with NLL^+z^ (dashed 
red lines). Pure soft contributions beyond 0{as) amount to 5 — 60% of the fixed-order 
NLO result, depending on the mass and process considered, whereas pure Coulomb effects 
and interference of soft and Coulomb corrections can amount to up to 60%. An exception 
to this is the squark-squark production process, where the effect of Coulomb corrections is 
small. This particular behaviour originates from cancellations between the cross sections 
for same-flavour squark production, where the repulsive colour-sextet channel is numeri- 
cally dominant and gives rise to negative 0(a^ In^ /3//3) corrections, and different-flavour 
squark production, where the corresponding term is positive, due to the dominance of the 
attractive colour-triplet channel. 

For squark-antisquark, squark-gluino and gluino-gluino production, a significant portion 
of the total Coulomb and soft-Coulomb corrections originates from bound-state effects 
below threshold. These correspond to the difference between the NLL and NLLno bs (dot- 
dashed purple) curves in the plots. For squark-antisquark and squark-gluino production 
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Figure 7: NLL i^-factor for squark-antisquark (top-left), squark-squark (top-right), squark- 
gluino (centre-left) and gluino-gluino (centre-right) production at LHC with ^/s = 7 TeV, 
and for the sum of the four processes (bottom). The plots show K^ll as a function of 
M for different NLL approximations: NLL (solid blue), NLLno bs (dot-dashed purple) and 
NLL^ (dashed red). See the text for explanation. 
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Figure 8: NLL i^-factor for squark-antisquark (top-left), squark-squark (top-right), squark- 
gluino (centre-left) and gluino-gluino (centre-right) production at LHC with = 14 TeV, 
and for the sum of the four processes (bottom). The plots show -R'nll as a function of 
M for different NLL approximations: NLL (solid blue), NLLno bs (dot-dashed purple) and 
NLL^ (dashed red). See the text for explanation. 
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Figure 9: NLL i^-factor for the total SUSY production rate at LHC with ^/s = 7TeV as a 
function of the gluino mass nig and average squark mass rriq. The dashed hne corresponds 
to the most recent exclusion limit presented in [2]. 




600 800 1000 1200 1400 1600 1800 2000 1.0 1.5 2.0 



Figure 10: Ratio of the NLL production-cross sections for the processes (12. 7p to the total 
NLL rate of coloured sparticle production ctsusy for the LHC with ^/s = 7 TeV. Left: Mass 
dependence for a fixed mass-ratio rriq = nig = M , Right: Dependence on the ratio nig/niq 
for a fixed average mass {niq + nig)/ 2 = 1.2 TeV. 
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bound-state corrections amount to 2 — 10% of the fixed-order NLO cross section, whereas 
for gluino-gluino production they can be as large as 30%. 

Figure [9] shows the NLL i^-factor for the total SUSY production rate at the 7 TeV LHC 
as a contour plot in the {rrig, mg)-plane. The r-dependence of the total resummed cross 
section arises from an interplay of the r-dependence of the single-process cross sections 
and of the relative dominance of the four subprocesses for a given r. The largest K- 
factor is obtained for rriq ~ 2 TeV and rrig ~ 1.4 TeV, with corrections of 50% to the NLO 
cross section. The plot shows also the recent exclusion limit published by the ATLAS 
collaboration in [2j assuming a simplified model of a massless neutralino, a gluino octet 
and degenerate squarks of the first two generations, while all the other supersymmetric 
particles, including stops and sbottoms, are decoupled by giving them a mass of 5 TeV. 
The limits are therefore not directly comparable to our results which treat the sbottom 
as degenerate with the light-flavour squarks, but are shown here as an indication of the 
current LHC reach. We do not attempt to estimate how resummation would affect the 
determination of this limit. However, one can observe that in the large squark-mass region 
the exclusion limit crosses regions with a i^'-factor bigger than 1.3, where resummation 
effects on the limit extraction might be relevant. 

Given the large effect of resummation, especially for squark-gluino and gluino-gluino 
production, it is interesting to study how the relative contribution of the four production 
processes to the total SUSY production rate is modified by the inclusion of NLL corrections. 
This is shown in Figure [TOl The qualitative behaviour of the relative contribution of the 
four different processes is very similar to the LO result (Figure [1]). However at large 
masses one can notice an enhancement of the squark-gluino production rate compared to 
the squark-squark channel (left plot), as one would expect from the larger NLL if-factor 
for the first processes. For a fixed average squark and gluino mass of 1.2 TeV (right plot) 
the relative ratios are basically unchanged for moderate values of r = rrig/mq, though one 
observes a significant enhancement of the squark-antisquark cross section for large gluino 
masses (r = 2). 

As a result of including the threshold-enhanced higher-order corrections, one expects 
that the uncertainty due to missing perturbative corrections is reduced compared to the 
NLO results. While for NLO the theoretical uncertainty arises from scale-variation only, 
the total theoretical error of the NLL results is obtained by adding scale and resummation 
uncertainties in quadrature, as defined in Section 12.41 The uncertainty bands for NLO, 
NLL^+h and NLL approximations are shown in Figure [TT]for the LHC with ^/s = 7 TeV and 
in Figure [12] for the LHC with ^/s = 14 TeV . In all plots the cross sections are normalized 
to unity at the central values of the scales and other input parameters. It is evident that 
the combined resummation of soft and Coulomb effects (NLL, solid blue) generally leads 
to a significant reduction of theoretical uncertainties compared to the NLO result (dotted 
black), especially in squark-antisquark and squark-squark production, where the error is 
reduced by a factor 2 or more in the large-mass region. The behaviour of NLL^+ft (dashed 
red) is more process-dependent, with basically no uncertainty reduction compared to the 
fixed-order NLO result for squark-antisquark production, and moderate effects for squark- 
gluino and gluino-gluino production. For squark-squark production (and, as a consequence 
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Figure 11: Total theoretical uncertainty of the NLO approximation (dotted black), full 
NLL resummed result (solid blue) and NLLs+/j (dashed red) at the LHC with y/s = 7 TeV. 
All cross sections are normalized to one at the central value of the scales. 



of the dominance of squark-squark production, for the total SUSY production rate) the 
uncertainties of NLL^+zi and NLL are very similar, due to the smallness of Coulomb effects 
in this particular production channel. The large reduction of the scale dependence for 
squark-antisquark production by soft-Coulomb interference effects is consistent with recent 
NNLL studies in this channel [22] that include the first Coulomb correction. 
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Figure 12: Total theoretical uncertainty of the NLO approximation (dotted black), full 
NLL resummed result (solid blue) and NLL^+z^ (dashed red) at the LHC with ^/s = 14 
TeV. All cross sections are normalized to one at the central value of the scales. 



3.2 Benchmark points for SUSY searches at LHC 

In addition to the grid files provided with the arXiv submission of this paper, we here 
present numerical predictions for some benchmark points at the LHC with 7 TeV centre- 
of-mass energy, in order to illustrate the effect of our NLL results on the production cross 
sections. We employ the sets of benchmark points defined in [31], that are compatible with 
recent LHC bounds and other data such as 6 — s'j, but not necessarily with constraints 
from the anomalous magnetic moment of the muon, or from the dark matter relic abun- 
dance. We consider the seven lines in the constrained MSSM (CMSSM) parameter space 
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defined in [31] and one line for tlie minimal gauge mediated SUSY breaking (mGMSB) 
scenario. For eacli line we selected one benchmark point expected to be relevant for 5 fb~^ 
of data and a second point relevant for 10 — 30 fb~ , where we naively extrapolate the reach 
of the 1 fb~^ LHC data [2ji3j that exclude CMSSM benchmark points and simplified models 
with total SUSY production cross sections of the order of ususy ^ 0.04 pb for rriq < rrig and 
c"susY ^ 0.1 pb for rrig > rrig. The mGMSB scenario we have selected has a quasi-stable 
neutralino as next-to-lightest SUSY particle (NLSP), so a similar reach as for CMSSM-type 
scenarios are expected. Since only the squark and gluino masses are relevant for the pro- 
duction cross sections, we have chosen points with a reasonable spread of masses and mass 
ratios, covering the range of average sparticle masses, 1.3 TeV < {rrig + mq)/2 < 1.5 TeV, 
and the mass ratios 0.75 ^ ^ ^ 1.12. This mass range is also compatible with the esti- 
mated discovery reach [50j of rrig ~ 1.3 TeV (1.6 TeV) for nig ~ rrig at 5 fb^^ (30 fb~^) 
and rrig ~ 0.8 TeV (1 TeV) for rriq ^ rrig in a CMSSM scenario with tan /3 = 0.45, Aq = 0. 
The remaining families of benchmark scenarios introduced in [3l] tend to have very similar 
mass ratios as our selected points. Therefore the relative contributions of the different 
production channels and the effect of the higher-order QCD corrections will be similar, 
although the decay chains and the resulting collider signatures can be very different. For 
some scenarios lighter squarks and gluinos than the ones considered here might still be 
allowed, for instance in GMSB with a stau NLSP. Predictions for such scenarios can be 
obtained by an interpolation of the grid files provided with the arXiv submission of this 
paper. 

The SUSY breaking parameters and the resulting mass spectrum of the coloured SUSY 
particles for the selected points is shown in Tables O |3] and H] together with our best NLL 
predictions for the total production cross section of light-flavour squark and gluino pro- 
duction (including simultaneous soft-gluon and Coulomb resummation as well as bound- 
state effects). Here rrig denotes the average mass of all squarks except the stops, fol- 
lowing the setup of [22l|5T]. The low-scale mass parameters have been generated using 
SUSY-HIT [52] employing SuSpect2 . 41 [53j with the standard model input rrit = 172.5 GeV, 
tts('^z) = 0.1172, rrib{rnb) = 4.25 GeVj^ The cross sections for the separate squark and 
gluino production processes are shown in Tables |5] and El The stops are always heavier then 
m^^ > 750 GeV for the considered benchmark points so direct stop-antistop production 
will be out of reach at the LHC with ^/s = 7 TeV (for some of the benchmark points, it 
might be possible to discover them in the gluino-decay products). Therefore we will give 
results for stop-antistop production separately in Section 13. 3[ 

For comparison, NLO results obtained using PROSPINO are also shown in the tables. Our 
setup agrees with the one in [T6|[5T] and the NLO results agree at the expected one-percent 
level or better with the results obtained from an interpolated grid using NLL-f ast [51]. For 
the NLO results the scale uncertainty is estimated by varying the factorisation scale in the 
interval M/2 < fif < M, with the renormalization scale set equal to the factorization scale. 

^ Note that in the mass-spectra quoted in [34] the masses are rounded to 5 GeV accuracy and only 
the hght-flavour squark masses are averaged. Furthermore version 2.3 of SuSpect has been used for the 
CMSSM benchmark points and SOFTSUSY for the mGMSB benchmark points. 
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Table 2: CMSSM benchmark points and total inclusive SUSY production cross sections 
for tan/3 = 10, Aq = 0. All masses are given in GeV. 
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Table 3: CMSSM benchmark points and total inclusive SUSY production cross sections 
for tan/3 = 40, Aq = —500 GeV. All masses are given in GeV. 



The tables also include the relative PDF uncertainties of the NLO cross sections. Based on 
experience with top-pair production [30], we expect that these agree with the relative PDF 
uncertainties of the NLL results at the relevant accuracy. Note that, due to correlations, 
the PDF uncertainty of the total SUSY production cross section is not equal to the sum 
of the uncertainties of the individual channels [I6] . For the NLL results we quote the total 
uncertainty including scale and resummation uncertainties, as discussed in Section 12.41 
In this case, the uncertainty of the total SUSY production cross section was obtained by 
neglecting correlations and adding the uncertainties of the individual production channels 
linearly, which gives a conservative upper bound. 

As can be seen from Tables Ej [3] and HI for the benchmark points with m^ < m^ the 
NLL correction to the inclusive squark and gluino production cross section is in the 15%- 
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Table 4: niGMSB benchmark points and total inclusive SUSY production cross sections 
for tan/3 = 15, iVmess = 1, -^mess = 10^ GeV. All masses and scales in GeV. 



range, which is consistent with Figure M For these points the theoretical uncertainty is 
reduced from ±14-15% at NLO to ±9-10% at NLL, in agreement with the behaviour seen 
in Figure [TTl The PDF uncertainty is also small for the benchmark points with rriq < rrig, 
since only the relatively precisely known quark PDFs are relevant at leading order for the 
dominant squark-squark production channel. For scenarios where the gluinos are lighter 
than the squarks, the NLL corrections grow to 20-30% and the theoretical uncertainty at 
NLL is at the 11-12% level. These features can be understood by considering the relative 
contributions of the different production processes in Figure [10] and Tables [5] and [61 For 
rrig < rrig, squark-pair production with moderate NLL corrections in the 6-10% range is 
the dominant production channel, with a non-negligible contribution from squark-gluino 
production with larger NLL corrections of the order of 30%. For rrig < rriq the roles 
are reversed, resulting in larger NLL corrections to the total rate and a slightly larger 
theoretical uncertainty, as can be seen from the uncertainties of the different production 
channels in Figure [TTJ Due to the large uncertainties in the current gluon PDF sets for 
large x, the PDF error for these points rises to the 10%-level. In gluino-pair production, 
especially for large gluino masses, the PDF error can even become ~ 30%. At a mass ratio 
^ < 0.75, as realized for the 40.3.1 and 40.3.5 points, gluino-pair production overtakes 
squark-pair production as the second-most important channel. For these points, the gluinos 
are relatively light, and the NLL corrections to gluino pair production are at most 45%. 
For benchmark points with heavier gluinos, the NLL corrections can grow to 70-80%, but 
the gluino-pair production rate is negligible for these points. For the moderate mass ratios 
considered here, squark-antisquark production is always much suppressed, but would have 
the second-largest cross section for — — )■ 2 . 

3.3 Results for stop-antistop production 

In this section we present results for resummation for stop-antistop production PP — > tjtj. 
At LO the cross section depends only on the stop mass, while at NLO it also presents 
a (much smaller) dependence on the mixing angle and the squark and gluino masses. 
For definiteness, we follow [I6] and fix all the parameters, except for the stop mass, to 
the values corresponding to the CMSSM benchmark point 40.2.4 with SUSY breaking 
parameters mo = 700 GeV, mi/2 = 600 GeV, Aq = -500 GeV and tan/3 = 40 [31]. Using 
SUSY-HIT [52] we obtain the relevant input parameters Trig = 1386 GeV, rrig = 1358 and 
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Table 6: NLL results for CMSSM benchmark points for tan/3 = 40, Aq = —500 and 
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Figure 13: -R^nll for stop-antistop production at the LHC at y/s = 7TeV (left) and y/s = 
14TeV (right) for different NLL approximations: NLL (sohd blue), NLLno bs (dot-dashed 
purple) and NLL^ (dashed red). 
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Figure 14: Total theoretical uncertainty of the NLO approximation (dotted black), full 
NLL resummed result (solid blue) and NLL^+z^ (dashed red) for titi at the LHC with 
^/s = 7 TeV (left) and ^/s = 14 TeV (right). All cross sections are normalized to one at 
the central value of the scales. 

cos Of: = 0.39. The numerical dependence on these parameters has been studied in [15] and 
found to be negligible. As discused in Section [273] we will limit ourselves to results for the 
lighter mass-eigenstate ti. 

Numerical results for the lightest stop mass eigenstates for masses in the 100-400 GeV 
range are presented in table [71 The NLL i^'-factor is plotted as a function of the mass in 
Figure [131 For the stop-mass range considered here the resummation effects are in the 15- 
25% range at the LHC with ^/s = 7 TeV and at the 12-20%-level at a centre-of-mass energy 
of 14 TeV. The NLL-corrections are therefore moderate for the mass-ranges accessible at the 
LHC, but larger than for light-flavour squark-antisquark production for the same masses 
(c.f. figure [7|). Contrary to the latter, for stops the most substantial contribution is given 
by pure soft resummation, with Coulomb effects, including bound-state corrections, in the 
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a{pp ^ titi)(pb), = 7 TeV 




NLO 


NLL 


ApDF 


-f^NLL 


100 


4.18+°!^ X 102 


4.77^°:«e^ X 102 


+2.0% 
-2.4% 


1.14 


200 


1.28+°;}^ X 10^ 


1.47^°:25 X 101 


+3.1% 
-4.1% 


1.15 


300 


1 90+0.17 


1 40+0.24 


+4.4% 
-5.3% 


1.17 


400 


2.12t'oil X 10-1 


2.50toil X 10-1 


+5.7% 
-6.4% 


1.18 



Table 7: Stop-antistop production cross sections. All masses are given in GeV. 

5% range of the NLO result. As the partonic cross sections for the gluon-fusion channel are 
identical for the light-flavour and top squarks, the different behaviour can be attributed 
to the quark-antiquark induced channel. For the light-flavour squarks it is dominated 
by S-wave colour- singlet production with a large, attractive Coulomb potential. For stop 
production the colour-singlet channel is absent, and only the P-wave colour-octet channel 
with a smaller repulsive Coulomb potential contributes. The increased relative size of the 
soft corrections compared to the light-flavour squarks follows from the P-wave suppression 
of the quark-antiquark channel and the resulting dominance of the gluon-fusion channel 
with larger soft corrections due to the colour factors Ca = Nc- 

The total theoretical uncertainty of the fixed-order NLO result and the resummed cross 
section, with and without Coulomb resummation, is compared in Figure [TH The width of 
the uncertainty band for NLL and NLL^+zi is similar, consistent with the observed small 
effect of Coulomb resummation. One can also notice that the resummed results shows 
almost no uncertainty reduction compared to the NLO result, except for the high-mass 
range. 

3.4 Comparison with Mellin-space results 

Results for NLL resummation of soft logarithms for squark and gluino production have been 
presented earlier in p^HT6] . These works adopt the so-called Mellin-space formalism, in 
which the threshold logarithms are exponentiated in Mellin-moment space, where singular 
terms appear as logarithms of the Mellin moment N, and the resummed cross section is 
then numerically inverted back to momentum space. In this section we compare these 
earlier predictions to the momentum-space formalism adopted here, using the numerical 
code NLL-f ast [SlJ to compute the Mellin-space resummed cross sections. Since NLL-f ast 
provides results for soft-resummation only (i.e. no Coulomb effects beyond 0{as) are 
included), for the comparison we introduce two additional NLL implementations: 

• NLL^: this implementation includes NLL resummation of soft logarithms, but no 
Coulomb or hard effects beyond 0{as)- This is achieved using Eqs. f l2.19p and ( I2.24p 
and setting /i^ = /i/. For the soft scale we adopt the running scale given in Eqs. 
(I228D, (^M- 
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• NLLs, fixed: as above, but with the running scale replaced by the fixed soft scale /i^ 
determined from Eq. fl2.27p . 

The i^-factor for NLL^ (solid blue), NLL,, gxed (dashed red) and the Mellin-space result 
(green dots) from NLL-f ast are shown in Figure [151 For clarity, we do not show the 
resummation uncertainties of the NLL^ and NLL,, gxed results that are significant, especially 
for smaller masses, as can be anticipated from the full NLL-results including Coulomb 
resummation in Figure El Even though the difference between NLL^ and NLLMeiim can 
be sizable, the envelope of the three curves in the plots in Figure [151 is well within the 
theoretical uncertainties of the predictions, and is thus correctly accounted for by our 
estimate of the intrinsic ambiguities of NLL resummation. 

The comparison of our default running-scale implementation NLL^ with NLLMeiUn shows 
different features for the various production channels. For squark-antisquark and gluino- 
pair production the agreement is overall reasonable, although the behaviour as a function 
of the mass is different in both cases, leading to a better agreement at large masses in 
the former case and at smaller masses in the latter. For squark-gluino production, the 
agreement also improves at larger masses. In contrast, for the squark-squark production 
channel, there is a constant shift, and the differences are sizable over the whole mass 
range. Since the NLL corrections are small for this production channel, the discrepancies 
in the cross-section predictions are however less important. For all processes, the fixed- 
scale momentum-space results NLLs,fixed are closer to the Mellin-space results, with similar 
magnitude and mass-dependence of the corrections. Since NLL^ and NLLMeiUn both resum 
threshold logarithms (In /3 and In respectively) appearing in the partonic cross sections, 
while NLLs^fixed resums logarithms in the hadronic cross section, this behaviour is some- 
what counter-intuitive and deserves further studies in the future|3 The spread of the three 
predictions for stop production is comparable to the one observed for the other SUSY pro- 
cesses in the same mass range. In this case, however, the two momentum-space predictions 
with a fixed and running soft scale show a better agreement with each other. 

It is interesting to note that the NNLL results for squark-antisquark production ob- 
tained in the A^-space approach in [22] are very similar to our best prediction for this 
channel in Figures [3 and [H Part of this agreement can be attributed to the fact that 
Ref. p2] includes the interference of the first Coulomb correction with higher-order soft 
corrections, which give a dominant contribution to our full NLL predictions (i.e. to the dif- 
ference between the red and blue curves in Figures [71 and [S]) , and to the difference between 
the NLL and NNLL results in [22] ■ However, our prediction includes higher-order Coulomb 
corrections and bound-state effects not included in the results of [22], while their results in- 
clude one-loop hard corrections and NNLL soft corrections not included in ours. Therefore 
the good agreement is to some extent fortunate, and cannot necessarily be expected for the 
other processes where the Coulomb corrections are of less relative importance (in particu- 
lar for squark-squark production), and where an NNLL analysis remains to be performed. 

^An analytic comparison of momentum-space resummation with fixed soft scale to A^-space resumma- 
tion has appeared recently for Drell-Yan production [54 , but no such investigation for the running scale 
with a lower cutoff has been performed yet. 
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Figure 15: Comparison of soft resummation in the moment-space formalism adopted in 
this work and in the conventional Mellin-space approach for squark-antisquark (top- left), 
squark-squark (top- right), squark-gluino (centre-left), gluino-gluino (centre-right), stop- 
antistop (bottom-left) production and the inclusive gluino and light-flavour squark cross 
section (bottom-right) at LHC with ^/s = 7 TeV. The plots show the i^'-factor for our 
default running-scale implementation (NLL^, solid blue), for the fixed-scale implementation 
(NLLs^fixed, dashed red) and for the Mellin-space result (NLLMciim, green dots) obtained 
with NLL-f ast. See the text for more details. 
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Furthermore, a combined NNLL soft/Coulomb resummation as performed for top-quark 
pair production in [30] would include higher-order Coulomb effect included neither in our 
present predictions nor in [22]. 

4 Conclusions and outlook 

We have performed a combined resummation of soft-gluon and Coulomb effects for all 
squark and gluino pair-production channels at the LHC, including bound-state effects, 
based on the formalism derived in [28l|29]. We have also extended the derivation of the 
factorization and resummation formalism to the case of stop-antistop production with a 
quark- ant iquark intial state, that proceeds through a P-wave. The corrections from higher- 
order Coulomb and bound-state effects and their interference with soft corrections are not 
included in previous predictions and can be sizeable, in particular for gluino-pair produc- 
tion, where they are as large as the soft-gluon corrections alone, and for squark-antisquark 
production where the effect is even bigger. For benchmark scenarios with moderate squark- 
gluino mass splitting the effect on the total inclusive squark and gluino production cross 
section is less pronounced but still relevant, with the total NLL corrections of the order 
of 15 — 30% of the NLO cross section. Therefore these predictions should be taken into 
account in the analysis of the upcoming LHC results expected this year. To facilitate the 
application of our results we provided numerical predictions for some of the benchmark 
points defined in [31], and include grid files with our results for squark and gluino masses 
in the 200-2000 GeV range at the 7 TeV LHC with the arXiv submission of this paper. 

Our results for soft-gluon resummation alone, obtained in the momentum-space resum- 
mation approach [^TH55] and with the scale-setting procedure introduced in [3U], agree 
within our estimate of resummation ambiguities with results obtained in the Mellin-space 
formalism [THUS], although a more detailed study of the relation between the approaches 
and the different scale choices within the momentum space framework would be desirable. 

In our analysis, the squarks and gluinos have been treated as stable, but contributions 
to the cross section from below the nominal production threshold have been included 
through bound-state corrections. This is expected to take the effect of small, but finite, 
widths to some extent into account. A more refined analysis is possible in our framework 
using a complex energy in the argument of the potential function in the factorization 
formula (12.131) . as done in recent studies of the invariant mass spectrum of gluino-pair 
and squark-gluino production [231425] . The investigation of the combined effect with soft 
resummation, as well as the extension to NNLL accuracy following Ref. [5U], is left for 
future work. 
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A NLL resummation for stop-antistop production 



In this appendix we derive the factorization formula fl2.13p for stop-antistop production 
from quark- ant iquark annihilation that proceeds in a P-wave state. While the aim is to 
establish resummation at NLL accuracy, the arguments suggest that the combined soft- 
Coulomb resummation can be performed at NNLL accuracy as well. Our reasoning applies 
to partonic subprocesses with a leading P-wave contribution, as relevant for stop-antistop 
production, but not to P-wave contributions to S-wave dominated subprocesses. 

Following [29j the scattering process qq ^ tt + X is described in terms of an effective 
field theory using soft-collinear effective theory (SCET) for the initial state quarks 

and potential non-relativistic QCD (pNRQCD) (SnHHS] for the stop and antistop. In this 
framework, the scattering amplitude in the full MSSM is expressed in terms of expecta- 
tion values of EFT operators multiplied by short-distance coefficients C, that contain the 
dependence on hard modes not present in the EFT (e.g. gluinos, hard off-shell gluons): 



A{qq-^itX) = J2 (ttX|CW'^"(/i)|gg) 



EFT ■ 



(A.l) 



The series in i accounts for the threshold expansion in powers of /3, while Ra = 1,8 
are the irreducible colour representations of the stop-antistop pair. Due to the threshold 
kinematics, (fci -|- A;2)^ ~ 4M^, no coUinear modes can appear in the final state Xg, that 
only contains soft modes. Both the short-distance coefficients C^^^'^" and the operators 
can carry open Lorentz and spin indices, that have been suppressed. For NLL 
resummation, only the colour-octet state is relevant |fl For the case of P-wave production, 
the leading i = term in the expansion of the scattering amplitude ( lA.ll) vanishes and the 
first non-vanishing term occurs for i = 1, i.e. it is suppressed by 0{P) compared to the 
S'-wave case. The P-wave production operator for a colour-octet state reads 



V2 



2 ^' 



X 



(A.2) 



Here the field ip"^ (x) is a non-relativistic field that creates the stop (antistop) while 
(^c) is a collinear (anti-collinear) field that destroys the initial-state quark with momentum 
~ V^/2n (antiquark with momentum ~ y/I/2n), with light-cone vectors n'^ = n'^ = 0, 
n ■ n = 2. The Wc (Wc) are collinear (anti-collinear) Wilson lines summing up collinear 
gluon emission to all orders [55l[56]. Dt denotes the projection of the covariant derivative 
D/^i = {d^ — igsT"'A1j^{x)) orthogonal to the timelike velocity of the stop-antistop pair. 



^ At NLO a finite, non-logaritlimic colour-singlet contribution is generated by real-gluon emission [46] 
that is suppressed by as/?^ compared to the leading P-wave contribution, and therefore not even relevant 
at NNLL. Since there is no leading S-wave contribution to qq ^ tt, the situation is simpler than for quarko- 
nium production, where singularities cancel among colour octet S-wave and singlet P-wave channels |63j . 
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At this stage, the operator (]A.2|) is written in NRQCD and contains the full gluon field 
A. This does not yet incorporate a systematic expansion in /3. For this purpose we adopt 
potential non-relativistic QCD (pNRQCD) [59H62], where only soft gluons with momenta 
scaling as 

q ~ (A.3) 
and the non-relativistic fields with so-called potential momenta scaling as 

ko ~ , Ml3 (A.4) 

are retained in the EFT, while potential gluons with momenta scaling like f lA.4p and modes 
with momenta of order M/3 are integrated out. After this step, the covariant derivative 
in ( 1A.2|) contains only the soft gluon field Ag with the same scaling as the soft momentum. 



A^ ~ (A.5) 

Therefore the term QsAs fj^j in the covariant derivative is suppressed by a factor QsP com- 
pared to the derivative term d^j ~ (3 and the production operator (]A.2p at leading power 
in pNRQCD is simply given by 



(/.). (A.6) 



The potential, collinear and anti-coUinear fields interact only through the exchange 
of soft gluons. These interactions can be removed from the leading effective-theory La- 
grangians by a decoupling transformation of the collinear [56|[M] and potential [2^ fields 
involving soft Wilson lines of the schematic form Wl^c = 5'„PVi''''^^c°^ and x = SwX^^\ and 
similar transformations for the conjugate fields. For definitions of the Wilson lines S'„ and 
5*^ see [25]. Since the fields with superscript (0) do not interact with each other in the 
leading effective Lagrangians, the scattering amplitude assumes a factorized form: 

Aiqq ^ klX) = mW^'^^6%Aq) (0| feW^.),, |g) 



Note that the factorization ( 1A.7|) holds only at leading power in the /3-expansion, since the 



field redefinition does not remove the interaction with the spatial components of the soft 
gluons in the P-wave production operators, or in the subleading effective Lagrangians. We 
will argue below that these corrections are not relevant at NLL accuracy for stop-antistop 
production from quark- ant iquark annihilation. 

Using the representation of the scattering amplitude flA.7p and following the steps 



in the derivation of the factorization formula discussed in detail in [29] one derives the 
factorization formula 



[s, /i) = J du J[{E - |) W\uj, /i) . (A.8) 
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Here the potential function for P-wave production in a colour octet state is given by 

J[iq) = \j d'ze'"-^ (0|[x(°)^S?T"^(°)](^)[^(°)^T'^K^X^°^](0)|0) , (A.9) 

while the colour-octet soft function (defined in terms of the squared matrix element of the 
soft Wilson lines in (IA.7I) ) is identical to that appearing for S'-wave production [28j. The 
factorization formula resembles that for S'-wave production, up to the replacement of the 
potential function by the appropriate P-wave expression. 

We now show that no corrections to flA.Sp appear at NNLL. As far as soft-gluon effects 
are concerned, we have argued that the corrections to the production operator flA.6p from 
the subleading A^^T-terms are of order QsP. Since the potential function for S-P-wave 
interference vanishes, there can be no contributions to the total cross section where one Ay- 
operator interferes with a leading operator (]A.6|) . The squared Ay-term gives a correction 
to the cross section of the ordeio Aa ~ a^as/S^jl, ln/3}. These corrections are at most of 
order N^LL in the combined counting as/P ~ 1 and a^ln/? ~ 1. 

Following [2^ one sees that insertions of higher-order 5o/i(- co//mear interactions into the 
matrix element (lA.ip are at least suppressed by and therefore are beyond NNLL. 

For P-wave production it has to be shown in addition that no contributions to the i = 1 
operators in (lA.ip are generated from partonic channels with leading S-wave contribu- 
tions. In the EFT language they would correspond to operators where one or both of the 
quark/antiquark fields are replaced by (anti)collinear gluon fields. A non- vanishing matrix 
element of such an operator in the \qq) intial state requires splittings of a collinear quark 
into a collinear gluon and a soft quark that are mediated by C(/3)-suppressed interactions 
in the SCET Lagrangian [58] (recall that no collinear final state particles appear at par- 
tonic threshold). Since there is no g^f- initiated production of stop-antistop pairs at leading 
order, two subleading splittings are required so these contributions to the amplitude are 
suppressed by ggf3 compared to the i = 1 term. The resulting real corrections to the total 
cross section are therefore even further suppressed compared to the subleading soft effects. 

Finally, subleading soft-potential effects due to chromoelectric x ■ Eg interactions [62] 
have been shown to be beyond NNLL in the S-wave case [29] since the potential function 
with a single insertion of an operator ~ x vanishes by rotational invariance. The same 
argument holds for the P-wave potential function!^ Furthermore, as mentioned before, due 
to the absence of a leading S-wave contribution to the quark- ant iquark channel, no mixing 
of S-wave and P-wave states by subleading interactions appears. Therefore, while we didn't 
perform an exhaustive study of effects beyond NLL, we do not encounter an obstruction 
for NNLL resummation for the P-wave contribution to stop-antistop production. Note, 
however, that Coulomb resummation at this accuracy requires the computation of the 
NLO P-wave Coulomb Green's function for scalar particles. 

®We have allowed for logarithmic corrections in the generalized potential function from contractions of 
the two gluon insertions. A simple estimate indicates, however, that these are in fact absent. 

^"^A possible interference of a single j4T-term in the production operator with a chromoelectric vertex 
might lead to a non- vanishing potential function but would be of the same order as/3^{l, ln/3} as the soft 
contributions considered above. 
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B Resummation functions 

The resummation functions appearing in f l2.18p are given by 



n / X exp 



(B.l) 



and T] = 2ar(/is, /i/). The labels r,r' denote the colour representation of the initial-state 
partons p,p'. At NLL the functions S, a( and or are given by 



Cr + Cr' 



(3o 



4:71 f 1 , 

— - — - 1 In r 



Po I 2po 



Inr, = ^-r— Inr, a^'^'i^i, fij) 



2/3o 



2/3o 



Inr, (B.2) 



where r = and 'jY = Y + Y' + Ih^s- The one-loop anomalous-dimension 

coefficients appearing here are given by 



{0),Ra _ cyf^ 



^R^. 1^'^' = -SC^ = -7^°)^'^ 7^°)' = -/3o = -7 
and the coefficients of the beta function are 



(0)0,8 



{3i = y - ^CaUj - 2CFnf. 



(B.3) 



(B.4) 



We also used the factor K = — Ca — ^TpUf appearing in the ratio of one- and 

two-loop cusp anomalous dimensions. The explicit values of the Casimir invariants for the 
SU (3) representations relevant for squark and gluino production are: 



4 10 
Ci = , C'a = - , ^6 = 



Cl0 = 6, C'i5 = y, C27 = 8. (B.5) 



C Determination of Hs and /3cut 



The soft scale /i^ used in the fixed-scale implementation, NLLfixed, is defined by Eq. fl2.27p . 
For equal squark and gluino masses, rrtg = rrig, the minimization procedure yields the 
following results: 



qq : 




= 127- 


170 GeV {y^ = 


7 TeV) 




= 146 - 


376 GeV (v^ = 


14 TeV) 


qq : 


f^s 


= 122 - 


135 GeV (v^ = 


7 TeV) 


fJ'S 


= 146 - 


312 GeV (v^ = 


14 TeV) 


M ■ 


f^s 


= 109 - 


141 GeV (v^ = 


7 TeV) 


fJ'S 


= 131 - 


310 GeV (v^ = 


14 TeV) 


99 ■ 


l^s 


= 106 - 


141 GeV (v^ = 


7 TeV) 


fJ'S 


= 127- 


308 GeV (v^ = 


14 TeV) .(C.l 
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LHC(7 TeV) 


LHC(14 TeV) 


qq 


a = 12.1328 
b = 159.828 


a = 11.5289 
b = 171.856 


qq 


a = 10.9627 
b = 265.414 


a = 10.9981 
6 = 278.062 


M 


a = 16.5459 
b = 235.227 


a = 15.3779 
6 = 267.763 


99 


a = 17.95 
b = 231.545 


a = 16.7483 
b = 264.702 


it 


a = 16.8468 
6 = 365.735 


a = 15.7282 
b = 481.045 



Table 8: Coefficients of the fit (IC.3I) for the squark-gluino production processes at centre- 
of-mass energies of 7 and 14 TeV. 



The scales for a centre-of-mass energy of 7 TeV refers to a mass range M = 500 — 2000 GeV, 
while for a 14 TeV LHC the mass interval M = 500 — 3000 GeV was considered. For stop- 
antistop production we obtain the scales 



tt 



fis = ^9- 138 GeV (^i = 7 TeV) /x. = 31 - 230 GeV (v^ = 14 TeV) (C.2) 



for the mass range = 100 — 1000 GeV for a centre-of-mass energy of 7 TeV and = 
100 — 1400 GeV at 14 TeV . In [33j it was suggested to fit the mass and energy dependence 
of the soft scales ( IC.lll by a function of the form 

Mjl-p) 
ya + bp 

with p = AJVP/ s. The coefficients a and b for the different processes are given in Table [HI 
and provide a fit to fIC.ip with accuracy better than 5% for the mass range considered. 

The procedure used to determine (3cut, which enters the definition of the soft scale used 
in our default NLL implementation, Eqs. ( 12.28^ and ( ]2.29p . was explained in detail in |30j . 
Following [30] we introduce eight different cross sections 



cut J 



^(/3cut - P) 



0"^ 



e{f3 - A 



cut J 



(C.4) 



defined using one of two possible matching prescriptions for the lower interval (y4< G 
{NLLi,NLL2}) and one of four possible approximations for the upper interval (i?> G 
{NLL2, NLOapp, NNLOa, NNLOb}). Here NLL2 denotes our default approximation, Eq. 
f l2.20p . while NLLi is the resummed result matched to the Born instead of the NLO cross 
section, 

(C.5) 



j.NLLi 
PP' 



-NLL/A\ 



NLL{0) , 



PP 
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Figure 16: Mass dependence of the parameter /3cut for the four SUSY production processes. 

NLOapp represents the sum of the full Born cross section and the approximated NLO 
corrections given in fl2.10p . while the two NNLO approximations contain in addition the 
0{a^) terms arising from the expansion of the NLL resummed result, including all of them 
(NNLOa) or only the subset which is completely determined at NLL (NNLOb). Pent is 
then determined such that the width of the envelope of the eight different implementations 
att{A^, 5>, /3cut) is minimal0 

The procedure described provides, for the default choice ks = 1 and equal squark and 
gluino masses, the following values for f3cut- 



qq : 


Pent 


= 0.47 


-0.36 (Vs = 


7TeV) 


Pcut 


= 0.48 


-0.38 (v^ = 


14TeV) 


qq : 


Pent 


= 0.30 


-0.24 (v^ = 


TTeV) 


Pcut 


= 0.31 


-0.25 (v^ = 


14TeV) 


M ■ 


Pcut 


= 0.40 


-0.34 (v^ = 


TTeV) 


Pcut 


= 0.40 


-0.36 (v^ = 


14TeV) 


99 ■ 


Pcut 


= 0.50 


-0.39 (v^ = 


TTeV) 


Pcut 


= 0.52 


-0.40 (v^ = 


14TeV) 



The numbers refer to the usual mass range of 500 — 2000 GeV at ^/s = 7 TeV and 500 — 
3000 GeV at a/s = 14 TeV, and the exact mass dependence is plotted in Figure UHl For 
stop-antistop production the result is 

ii: /3cut = 0.53 - 0.40 (v^ = 7 TeV) /3cut = 0.54 - 0.41 (^^ = 14 TeV), (C.7) 

again for the mass range mj = 100 — 1000 GeV at 7 TeV and = 100 — 1400 GeV at 
14 TeV. The theoretical ambiguities of the resummed result with the running scale are 
estimated as follows: i) the default implementation NLL2 is parametrized in terms of E 
instead of (3, ii) the NLL2 result is computed for the values kg = 1/2,2, recomputing /3cut 
anew for each choice and iii), the envelope of the eight cross sections flG.4p is taken while 

^^The values (|C.6p are also used for the NLL^+^i approximation (|2.19|) . instead of recalculating /3cut 
using that approximation in (|C.4p . Since the difference of NLL and NLLs+h is used in order to assess the 
effect of Coulomb resummation, we consider it more meaningful to use the same /3cut for both in order not 
to obscure the genuine Coulomb effects by different scale choices. 
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/Scut is varied by ±20%. The uncertainties from the three sources are added in quadrature, 
and iii) gives the dominant contribution, in the 3 — 13% range depending on the process, 
while i) and ii) are usually below 2%. 
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